Tracing the redshift evolution of Hubble parameter with gravitational-wave standard 

sirens 



O 

5-1 

< 

o 

u 

43 

9 L1, 
6 

03 



(N 
> 
O 
O 

o 
in 



Atsushi Nishizawa/'B Atsushi Taruya, 2,3 and Shun Saito 4 ' 5 

1 Yukawa Institute for Theoretical Physics, Kyoto University, Kyoto 606-8502, Japan 
2 Research Center for the Early Universe, Graduate School of Science, 
The University of Tokyo, Tokyo 113-0033, Japan 
3 Institute for the Physics and Mathematics of the Universe, 
The University of Tokyo, Kashiwa, Chiba 277-8568, Japan 
^Department of Physics, Graduate School of Science, 
The University of Tokyo, Tokyo 113-0033, Japan 
5 Department of Astronomy, 601 Campbell Hall, University of California Berkeley, California 94720, USA 

(Dated: April 28, 2011) 

Proposed space-based gravitational-wave detectors such as BBO and DECIGO can detect ~ 10 6 
neutron star (NS) binaries and determine the luminosity distance to the binaries with high precision. 
Combining the luminosity distance and electromagnetically-derived redshift, one would be able to 
probe cosmological expansion out to high redshift. In this paper, we show that the Hubble parameter 
as a function of redshift can be directly measured with monopole and dipole components of the 
luminosity distance on the sky. As a result, the measurement accuracies of the Hubble parameter 
in each redshift bin up to z — 1 are 3 — 14 %, 1.5 — 8 %, and 0.8 — 4 % for the observation time 1 yr, 
3yr, and 10 yr, respectively. 

PACS numbers: 



I. INTRODUCTION 

Future space-based gravitational wave detectors such 
as DECI-hcrtz Interferometer Gravitational-wave Obser- 
vatory (DECIGO) 0, II and Big-Bang Observer (BBO) 
Q (also see Q for updated information) are the most sen- 
sitive to a gravitational wave (GW) in 0.1 — 1 Hz band 
and will aim at detecting a cosmological GW background 
generated during the inflationary epoch, the mergers of 
intermediate-mass black holes, and a large number of 
neutron star (NS) binaries in an inspiraling phase. These 
GW sources enable us to measure the cosmological ex- 
pansion with unprecedented precision to investigate 
the population and formation history of compact binary 
objects, and to test alternative theories of gravity [1, H[. 

It has been known that the continuous GW signal from 
a compact-binary object provides a unique way to mea- 
sure the luminosity distance to the source with high pre- 
cision. Such binary sources are often referred to as stan- 
dard siren (analogous to the electromagnetic standard 
candle). With the redshift information determined by 
an electromagnetic follow-up observation, the standard 
siren can be an accurate tracer of the cosmic expansion 
Q. The feasibility of the standard siren relics on the 
determination of the redshift of each binary. The identi- 
fication of a host galaxy by follow-up observation arc thus 
crucial, indicating that a high-angular resolution is gen- 
erally required for GW detector. In the case of DECIGO 
and BBO, the detectors orbit the Sun with a period of one 
sidereal year, and constitute four clusters, each of which 
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consists of three spacecrafts exchanging laser beams with 
the others. The two of the four clusters are located at 
the same position to enhance correlation sensitivity to 
a gravitational wave background, and the other two are 
separated on the Earth orbit in order to enhance the an- 
gular resolution so that we can easily identify the host 
galaxy of each NS binary via the electromagnetic follow- 
up observations. Based on this setup, Cutler and Holz 
[J] have shown that cosmological parameters can be accu- 
rately measured by DECIGO and BBO with a precision 
of ~ 1%, assuming a large number of neutron star (NS) 
binaries, ~ 10 6 . 

In this paper, we show that the space-based GW detec- 
tors can also measure the Hubble parameter H(z) from 
the GW standard sirens (3f|. I n this method, the mea- 
sured quantity is independent of that in the usual method 
of a standard siren, in which an observable is the lumi- 
nosity distance as an integrated quantity of H~ 1 (z). The 
method to measure the Hubble parameter at each red- 
shift has been proposed by Bonvin et al. [1, [§j], who 
originally developed this idea in the observation of dis- 
tant type la supernovae. In general, a large number of 
samples is necessary for the accurate measurement of 
H(z), and they concluded that 10 5 - 10 6 supernovae are 
needed to achieve a few percent accuracy. In contrast to 
the supernovae observation which requires unrealistically 
large number of the samples and suffers from relatively 
large systematics, DECIGO and BBO are expected to 
detect a million of NS binaries with smaller systematic 
errors. Thus, the measurement of the Hubble parameter 
at high redshifts becomes even more feasible with stan- 
dard sirens. 

The Hubble parameter H (z) can be also measured by 
estimating the differential age of the oldest galaxies in 
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each redshift bin and using the baryon acoustic os- 
cillation (BAO) along the line-of-sight direction from the 
spectroscopic galaxy samples (e.g., [HKlH). ^ n this re- 
spect, the present method with GW observation is com- 
plementary, and useful for an independent cross check. 

This paper is organized as follows. In SccHH we briefly 
review the basic idea to measure the Hubble parameter 
via dipole of the luminosity distance g?l, originally pro- 
posed by Bonvin et al. H Q , and derive the basic equa- 
tions to estimate the accuracy of the Hubble parameter. 
In Sec. lIIIl the method is applied to the GW observation. 
We briefly describe the GW standard siren and calculate 
the measurement accuracy of the Hubble parameter. In 
Sec. IIV1 systematic errrors are discussed and compared 
with the uncertainty coming from instrumental noise in 
GW observation. Finally, Sec. [V] gives a brief sum- 
mary and discussions on the feasibility of our method. 
Throughout the paper, we adopt units c = G = I, and 
assume a flat universe. 



II. HUBBLE PARAMETER FROM THE DIPOLE 
ANISOTROPY OF LUMINOSITY DISTANCE 

Consider the luminosity distance to some astronomi- 
cal objects measured at redshift z and angular position 
n. In principle, the observations of many objects over 
the sky enables us to map out the angular distribution 
of luminosity distance, di(2,n), and no directional de- 
pendence appears if the observer is at cosmological rest- 
frame (i.e., CMB rest- frame) in a homogeneous universe. 
However, there certainly exist tiny anisotropics in cLl aris- 
ing from the matter inhomogeneities of the large-scale 
structure and/or the local motion of the observer [R|. As 
it has been shown in Ref. @ , the dominant component of 
anisotropics is the dipole induced by the peculiar velocity 
of the observer, and the contribution to the higher multi- 
poles coming from the weak gravitational lensing effect is 
basically small [36] . Then, we can expand the luminosity 
distance as 



d L (z,n) = df\z) + d^'(z) cos 



(i) 



— J dn d L (z,n), 



:i) (z) = JL J dn(n-e)d L (z,n), 



(1) 



(2) 



where we define cos 8 = n • e, and the quantity e is the 
unit vector directed toward the dipole. 

In the expression ((T|), the first term in the right-hand 
side is the direction-averaged luminosity distance, which 
is identified with the one defined in the homogeneous and 
isotropic universe: 



df(z) = (l + z) 



dz' 



(3) 



with the Hubble parameter given by 
H(z) = H {f2 m (l + zf + (1 — O ro )(l + z )^+^+ Wa ) 



x exp 



-3w a 



l + z 



1/2 



(4) 



Here we assumed the flat universe, and the dark energy 
equation-of-state parameter w(a) = P/p parametrized 
by w(a) = wq + w a (l — a). For later analysis of the 
error estimation, we adopt a fiducial set of cosmological 
parameters: Q m = 0.3, wq = —I, w a — 0, and Ho = 
72 km s" 1 Mpc" 1 . 

On the other hand, the second term on the right-hand- 
side of Eq. ([T]) arises from the Doppler effect due to the 
motion of the observer. To derive the expression for d^\ 
we approximate the propagation of GW or light follows 
the trajectory of a null geodesic. Writing the luminosity 
distance to an astronomical object as a function of the 
conformal time rj at which the source emits GW or light, 
the Doppler effect leads to [1, Q 



d L (r,,n)=d ( £\ V )[l-n-v ] , 



(5) 



where the vector vo indicates the peculiar velocity of 
a local observer, which is small enough relative to the 
light velocity. Note that the motion of the local ob- 
server also induces the Doppler effect in the redshift, 
z = z + Sz, where the unperturbed redshift z is defined 



as z = a (i]) — 1. To first order in vo, we have 
<Jz = -(n-v )(l + z). 



(0) 



Rewriting Eq. ([5]) in terms of the redshift z, then expand- 
ing d^ (z — Sz) up to the first order, and using Eqs. (|3]) 



and 



give 



d L (z,n) 



ddf{z) 
~^z~ ( 



df{z) 



[I - n • v ] 



4 0) (,) + 



H{z) 



(n ■ v ) 



(7) 



Comparing the above expression with Eqs.([T]) and ©, 
we arrive at 



d { l\z) 



|vol(l + *) S 
H(z) 



(8) 



with the direction of the dipole specifically chosen as 
e = vo/|vo|. In the expression ([8]), the magnitude of 
dipole anisotropy is inversely proportional to the Hub- 
ble parameter at the source redshift, because the per- 
turbed luminosity distance corresponds to the deriva- 
tive of Eq. ([3]) . Recalling the fact that the motion of 
the local observer also induces the same size of dipole 
anisotropy in the CMB and its amplitude is estimated as 
| v | = 369.1 ± 0.9km s" 1 [TJ]. Then the dipole anisotropy 
in the luminosity distance to high-z objects gives a direct 
measure of H(z). 
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Now, let us discuss the statistical error in the measure- 
ment of the Hubble parameter. We add the measurement 
error of the luminosity distance 6d,L(z,n) into Eq. ([T]): 
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d L (z, n) = df{z) + d^>(z)(n • e) + 5d L (z,n) . 

Then the measurement error of d^ is estimated from the 
definition of d^' , and is expressed as 

5d^\z) = -j^- J dn(e ■ n)SdL(z,n) . 

The variance of this becomes 



(i) 



Ad^(z) 



dn / dn' 
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x (e ■ n)(e • ri){5d L (z, n)5d L (z, n')) , 

(9) 

where the bracket represents the ensemble average over 
the sources given at z. Assuming that the distance errors 
measured from each GW source are statistically indepen- 
dent and isotropic, we have 

(Sd L (z,n)Sd L (z,n')) = 4tt \Adf\z)} 2 S 2 (n - n') 



where the quantity 



Adg>(*) 



is the variance of d 



(0) 



Substituting this into Eq. ©, we obtain 
2 9 r 



Ad£\z) 



47T 
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Ad[°\z) 
Ad*\z) " 



dn (e • n) , 



(10) 



From this and the expression ([8]) , the measurement error 
of the Hubble parameter at a given redshift z is related 
to the error of the direction-averaged luminosity distance, 
Ad]* (z), and is given by 



AH(z) 



Ad^(z) 



H(z) 4D ( z ) 



= V3 





-1 


~Adf{z) 






[ df{z) \ 



(11) 



Here we ignored the velocity error in CMB observation. 
In fact, the contribution is negligible as we will show in 
SecHV] The factor [d^ /d^}' 1 is typically very large 
in the present case with |vo| <C 1. In Fig. [TJ we plot 
the ratio of the dipole to the monopole in the luminosity 
distance. Figure implies that even a negligibly small error 

in the averaged distance d!f^ can produce a large scatter 
in H(z). 



10"- 



10 




FIG. 1: Ratio of (2) to df } (z), in which the CMB dipole of 
\vo\ = 369.1 ± 0.9 km/sec measured by Wilkinson Microwave 
Anisotropy Probe (WMAP) [3 is used. 



Note that the mean error on the Hubble parameter is 
reduced to AH/ y/N if we observe N independent sources 
at a given redshift bin. Thus, for a precision measure- 
ment of Hubble parameter, we need not only an accu- 
rate determination of averaged distance to each source, 
Ad^ /d^ <C 1 but also a large number of sources. In the 
case of the type la supernovae as standard candles, the 
averaged distance error is related to the intrinsic mag- 
nitude error Am as Ad^/d^ = (In 10 /5)Ato. Adopt- 
ing an optimistic value Am = 0.1, it yields the error 
Ad^ /d^ ~ 0.05, but to achieve a few percent accuracy 
in the Hubble parameter, we need unrealistically large 
number of samples of the order 10 6 . In other words, for 
a reasonable number of 10 4 samples (sec e.g. [l^, [la]), 
the systematics in the averaged distance should be re- 
duced to Am < 0.01, which seems very difficult from the 
empirical calibration. 



III. GW STANDARD SIRENS 



In this section, we consider the GW standard sirens 
as an alternative probe to measure the Hubble param- 
eter from the dipole anisotropy of luminosity distance. 
The advantage to use the standard sirens observed by 
space-based GW detectors is that the expected num- 
ber of sources (NS binaries) is much larger than that of 
the type la supernovae, and the sources are distributed 
deeply enough at higher redshifts. Moreover, the NS bi- 
nary would be a clean GW source, and with DECIGO or 
BBO, the luminosity distance d^z, n) can be accurately 
measured with less systematics. 
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A. Luminosity distance error 

Let us first estimate the distance error of the stan- 
dard siren, taking account of the instrumental noise of 
the GW detector. Possible systematic errors will be dis- 
cussed later. In GW experiments, a direct observable is 
the waveform of the GW signal, and comparing it with a 
theoretical template, we not only determine the system 
parameters of GW source but also extract the cosmolog- 
ical information. 

For a single binary system, the Fourier transform of the 
GW waveform is expressed as a function of frequency / 

0,(3, 



kf) 



A 



d L (z) 



7/6 e i*(/) 



(12) 



where is the luminosity distance, and the quantity 
M z = (1 + z)77 3//5 M t is the redshifted chirp mass with 
the total mass M t = mj + m 2 and the symmetric mass 
ratio r\ = mim2 / Mf . Here, the constant A is given by 
A = (V6 7r 2/3 ) _1 , which is multiplied by the factor \[kfh 
for a geometrical average over the inclination angle of a 
binary [37j . The function represents the frequency- 
dependent phase arising from the orbital evolution, and 
at the order of the restricted 1.5 post- Newtonian (PN) 
approximation, it is given by (l7l . Il8j 



*(/) 



27r/t c -&-- + — (ttMJ)- 5 ^ 
-167T77- 3 / 5 (7rM z /)l , (13) 



where t c and <f> c are the time and phase at coalescence, 
respectively. The first term in the bracket in Eq. (JT3J) 
corresponds to Newtonian-order dynamics and the other 
remaining terms represent the Post- Newtonian order cor- 
rections in powers of v ~ (nMzf) 1 / 3 . In principle, there 
additionally appears a phase correction due to cosmic 
expansion, and the Hubble par ameter H(z) can be also 
measured from this term (II. Il9j|. Although the inclusion 
of the phase correction slightly changes the size of the 
errors in binary parameters, it does not seriously affect 
the estimation of the luminosity distance d^. In addition, 
the sensitivity of the phase correction to the Hubble pa- 
rameter is rather small. Thus, we may safely ignore the 
phase correction due to cosmic expantion in the subse- 
quent analysis. 

In Eqs. (|12[) and (|13p , there are five unknown param- 
eters to be determined obscrvationally, i.e., M Zl 77, t c , <p c , 
and dh- Except for the luminosity distance, the four pa- 
rameters merely carry the information on the individual 
property of the binary system. For simplicity, we con- 
sider the equal-mass NS binaries with 1.4M , which lead 
to M z = 1.22(1 + z)M Q and 77 = 1/4, and set the other 
parameters to t c = and <j) c = 0. 



Since the GW observation can only determine the red- 
shifted chirp mass M z , the redshift of each binary has to 
be measured from an electromagnetic counterpart. Ac- 
cording to Cutler and Holz Q, the angular resolution of 
BBO is ~ 1 — lOOarcsec 2 , with which we can identify 
the host galaxy of the binary. We thus suppose that the 
redshift of any binary system is obtained from the electro- 
magnetic observations. Note that the Doppler effect by 
the local motion also affects the redshifted chirp mass, 
and the dipolc anisotropy might be measured through 
the spatial distribution of the observed chirp mass if 
the intrinsic scatter in the mass distribution of NS bina- 
ries is very small. The feasibility to measure the dipole 
anisotropy from the chirp mass might be interesting, but 
we need a more detailed study on the formation history 
of NS binaries, and we here simply ignore this effect in 
the parameter estimation. 

The fundamental basis to estimate the distance error 
for a single binary is the Fisher matrix formalism. The 
Fisher matrix for a single binary is given by [l7|, 

where d a denotes a derivative with respect to a parameter 
a ; M z , 77, t c , 4> c . and dh- The quantity represents the 
GW signal obtained from the z-th interferometer. Since 
two independent signals are obtained for each cluster [2ll ] , 
DECIGO has the eight interferometric signals in total, 
each of which is supposed to have an identical detector 
response and noise power spectrum P(f). The analytical 
fit of noise spectrum (38j is given by 



P(f) 



4.21 x 10" 



'(lHz) 



1.25 x 1CT 47 



+3.92 x IQ-^l-^-j Hz" 1 . 



In Fig. the noise spectrum of DECIGO is shown, to- 
gether with the evolutionary tracks of the NS binary lo- 
cated at three different redshifts, z = 0.1, 1, and 5. In 
each track, the symbols indicate the frequency at the 10, 
3 and 1 years before the time of binary coalescence (from 
left to right). In this respect, the lower cutoff of the fre- 
quency /min should be incorporated into the integration 
in Eq. (|L4")l. and is given by the function of observation 
time T bs as well as the redshift and mass: 



/min — 0.233 



\ 5/8 /-, \ 3/8 



M z 



Tobt 



Hz 



(15) 



Note that the coalescence frequency of the NS binary is 
typically ~ kHz, and thus the upper cutoff of the fre- 
quency naturally arises from the noise curve. For the 
computational purpose, we set / max = 100 Hz. 

Given the numerically evaluated Fisher matrix, the 
marginalized 1-sigma error of a parameter, A9 a is es- 
timated from the inverse Fisher matrix 



A6 a = ^{T^jaa. 



(16) 
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/[Hz] 

FIG. 2: Sky-averaged DECIGO noise curve. (Arm angle 60° 
is taken into account.) Diagonal lines represent frequency 
evolutions of a NS-NS binary at z — 5 (solid, red), 2 = 1 
(dotted, green), and z = 0.1 (dashed, blue). Diamonds on 
the lines from the right to the left denote the frequency of the 
binary lyr, 3yr, and 10 yr before the merger. 
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FIG. 3: Measurement accuracy of the luminosity distance 
with a single binary as a function of redshifts. The curves 
tagged (Jmst are those determined only by instrumental noise 
and with the observation time 1 yr (red, solid curve) , 3 yr 
(green, dotted curve), and 10 yr (blue, short-dashed curve), 
respectively. The lensing error and the peculiar velocity error 
are represented by magenta (long-dashed) and light blue (dot- 
dashed) curves. 



In Fig. |31 the resultant error of the luminosity distance 
for a single binary, <7i ns t, is plotted against a source red- 
shift, assuming the observation time lyr (solid curve), 
3yr (dotted curve), and 10 yr (dashed curve). The over- 
lap of these three curves indicates that <7i, ls t hardly de- 
pends on the observation time, because the observation 
time appears only through the cutoff frequency / m i n with 
the fractional power of 3/8 and the improvement of the 
precision is generally slow. Even for a single binary sys- 
tem, the precision of a few percent level is easily achiev- 
able for the distance measurement in the absence of sys- 
tematic errors, and this is also true for a rather high-z 
binary. 

B. Accuracy of Hubble parameter 

Given the uncertainty of the averaged luminosity dis- 
tance for each binary, the accuracy of the Hubble param- 
eter is estimated from Eq. (JTTJ) , and with the ensemble 
over the AN(z) independent binary systems in the vicin- 
ity of the rcdshift z, we can get an improved constraint 
on the Hubble parameter at each redshift bin. 

Here, to derive the measurement error of the Hubble 
parameter, we adopt the follo wing fitting form of the NS- 
NS merger rate given by Ref. |22| : 

{l + 2z (z<l) 
|(5- z) (1< z < 5) , 
(5 < z) 

(17) 

where the function r(z) is estimated based on the star 
formation history inferred from the UV luminosity [23| . 
The quantity no represents the merger rate at present. 



Though it is still uncertain, we adopt the most recent 
estimate, ho = 10~ 6 Mpc~ 3 yr -1 , as a reliable and confi- 
dent estimate based on extrapolations from the observed 
binary pulsars in our Galaxy [24| . Then, the number of 
NS binaries in the redshift interval [z — Az/2, z + Az/2] 
observed during T b s , AN(z), is given by [22j 

f z+Az/2 ./a 

AJV(z) = T obs / dV c (z') , (18) 

Jz-Az/2 i- + Z 

where dV c means the comoving volume element defined 
as dV c (z) = Airr 2 (z)dz / H (z) with the comoving radial 
distance r(z) = c?l(z)/(1 + z). 

In Fig. 21 observed redshift distribution of NS binaries 
AN{z) is plotted, assuming the 3 year observation and 
the redshift width Az = 0.1. The total number of NS bi- 
naries is ~ 10 6 , which is much larger than the expected 
number of type la supernovae. Note that the number 
of merger events increases with T b s , and thus the ac- 
curacy of the Hubble parameter is improved by a factor 

1/2 

-^obs • Combining this and the distance error in previ- 
ous subsection, Fig. [5] shows the expected errors for the 
Hubble parameter measured from the dipole anisotropy. 
The three different error bars in each rcdshift bin repre- 
sent the results from the 1-, 3-, and 10-year observations 
(from large to small sizes). The Figure implies that up 
to the redshift 2 = 1, the Hubble parameter can be ac- 
curately measured with a precision of 2 — 5%, 1 — 3%, 
and 0.7 — 1.5% for the observation time of 1, 3, and 10 
years, respectively. Even at z = 2, the Hubble parameter 
can be measured with a precision of 18%, 10%, and 6% 
for the observation time of 1, 3, and 10 years, respec- 
tively. This is quite impressive in the sense that a GW 
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FIG. 4: Number of NS-NS binaries (in the unit of 10 4 ) that 
would be observed by DECIGO in each redshift bin of Az — 
0.1 at a redshift z during 3yr observation. As is manifest 
from Eq. (|18[) . the number of the binaries scales linearly with 




FIG. 5: The Hubble parameter calculated with our fiducial 
cosmological parameters (solid curve) and la-error bars esti- 
mated in the cases that we use all binaries observed by DE- 
CIGO during the observation time, lyr (red), 3yr (green), 
and 10 yr (blue) (long observation time corresponds to the 
smaller error bar). 



standard siren has a nearly equal sensitivity to the Hub- 
ble parameter with other complementary methods such 
as BAO. Another noticeable point using the standard 
sirens is that we can trace the redshift evolution of Hub- 
ble parameter even at higher redshift z > 1. Although 
the number of high-z NS binaries is highly uncertain, the 
standard sirens would be potentially powerful to probe 
the early-time cosmic expansion, and should deserve fur- 
ther investigation. 



IV. SYSTEMATIC ERRORS 

So far, we have discussed the accuracy of Hubble pa- 
rameter taking only account of the distance error associ- 
ated with the instrumental noise. However, there are sev- 
eral effects which may systematically affect the measure- 
ment of dipole anisotropics in the luminosity distance, 
leading to increasing the error in the Hubble parameter. 
Among them, a dominant contribution may come from 
the gravitational lensing magnification induced by the 
matter inhomogeneities of large-scale structure along the 
line of sight (e.g., j25rl28| ). which systematically changes 
the luminosity distance to each binary system. Another 
important effect would be the peculiar velocity of the bi- 
nary along the line of sight, which randomly contributes 
to measurement error via Doppler effect. These system- 
atic errors to the averaged luminosity distance are sum- 
marized as 



Ac4 0) (z) 

l 4 0) (*) 



hist 



(19) 



where <7i ns t is the error associated with the GW experi- 
ment in Sec. IIII Al <r\ cns is the lensing error, and a pv is 
the peculiar-velocity error. 

There are several studies on the effect of lensing mag- 
nification, particularly focusing on the distance measure- 
ment from the type la supernovae. Holz and Linder 
[26l | estimated the lensing error on the distance measure- 
ment by using Monte Carlo simulation, and assuming the 
Gaussian form of lensing magnification probability, they 
derived a fitting formula for the systematic error. Later, 
the significance of non-Gaussian tail has been recognized 
0,[3O| , and it turned out that this effect reduces the lens- 
ing error by a factor of 1.5 - 2, compared to the Gaussian 
distribution. More recently, Hirata, Holz, and Cutler [29| 
adopted a log-normal distribution for the magnification 
probability and obtained the fitting formula for the (av- 
eraged) distance error: 



0"lcns(z) = 0.066 



l-(l + z) 



-0.25 



0.25 



1.8 



(20) 



In what follows, we adopt the lensing error in Eq. (|20p. 

As for the peculiar velocity error, the clustering of 
galaxies induced by the gravity leads to a coherent and/or 
virializcd random motion, which gives rise to the Doppler 
effect and affects the determination of cosmological red- 
shift via the spectroscopic measurement. In addition, 
binary barycentric motion itself in the host galaxy also 
leads to the Doppler effect, which causes random fluctu- 
ations in the luminosity distance. These two systematic 
effects can be of the same order and can be translated 
into the distance error as [3lf 



(1 + z) 2 



H{z)d L {z) 



^v,gal 
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FIG. 6: Measurement accuracy of the Hubble parameter with 
all observed binaries. We plot the measurement accuracies, 
including the lensing error, with the observation time T h s = 
1 yr, 3 yr, and 10 yr from the top to the bottom, respectively 
(solid curves), and those without the lensing error (dotted 
curves) . 



Here, cr Vjga i is the one-dimensional velocity dispersion of 
the galaxy. Taking account of the non-linear effect of 
gravity, it often set to ov^i = 300km s" 1 , mostly inde- 
pendent of the redshifts |32| • 

In Fig. [3j the distance errors for a single binary sys- 
tem from the lensing magnification, <7i en s, and peculiar 
velocity, cr pv , are overlaid, together with the uncertainty 
from the GW observations. It turns out that the lens- 
ing magnification could dominate the distance error, and 
exceeds the error from the GW experiments at z > 0.2. 
Thus, the lensing magnification could be potentially a 
main source of the distance error, and the accuracy of 
the Hubble parameter might be somewhat degraded. 

In Fig. |6j the size of the measurement error for Hub- 
ble parameter, AH(z)/H(z), is plotted as function of 
rcdshift, taking account of all the systematics and the 
binary distribution. Comparing with the case without 
lensing error (dotted curves), the resultant accuracy of 
Hubble parameter are degraded. Nevertheless, even in- 
cluding the lensing systematics, the Hubble parameter 
up to z < 1 can be accurately measured with the preci- 
sion 3-14 %, 1.5 - 8 %, and 0.8 - 4 %, for T obs = 1 yr, 
3yr, and 10 yr, respectively. Although the lensing effect 
is potentially crucial for the cosmological application of 
standard sirens, the technique to reduce the lensing ef- 
fect has been recently exploited [29|, IH, H3] , and the 
feasibility of the method has been discussed. With an 
improved technique developed near future, the lensing 
systematics would be removed, and one could approach 
the limit determined by the instrumental noise. 

We briefly comment on the uncertainty in the ampli- 
tude of dipole Vo . The current constraint on the motion 
of local observer comes from the CMB observation, and 



the estimated error of vo directly affects the accuracy 
of Hubble parameter in Eq. (|11|) . However, the current 
observation produces 

H(z) |v | 

Therefore, the systematics in the dipole from the CMB 
observation gives a tiny contribution, and can be ignored. 

Finally, note that the results obtained here rely on the 
fact that we can successfully identify the redshifts of all 
host galaxies. In practice, the identification would be 
difficult as increasing the redshift, because galaxies at 
high-z arc fainter and the time required for a spectro- 
scopic measurement is much longer than that for low-z 
galaxies. The resultant accuracy of the Hubble parame- 
ter me asurem ent would be degraded, being proportional 
to l/y/N(z), where N(z) is the number of galaxies in 
a bin at the redshift z. Thus, the measurement accu- 
racy depends on the capability of galaxy redshift survey 
available in the future. 

Another important issue is detector calibration, which 
would potentially affect the measurement accuracy of the 
luminosity distance. It is rather crucial not only for the 
Hubble parameter measurment, but also for the subtrac- 
tion of the neutron star binary foreground in order to 
achieve the detection of an inflationary gravitational- 
wave background as a primary science goal of DE- 
CIGO/BBO. Although this is beyond the scope of this 
paper, the issue should be considered seriously and ad- 
dressed in the future. 



V. SUMMARY 

In this paper, we have shown that the redshift evolu- 
tion of Hubble parameter can be measured by utilizing 
a large number of NS binaries observed by space-based 
GW detectors such as DECIGO and BBO. Although this 
method requires a large number of samples 10 5 -10 6 to 
measure H(z) with the accuracy of a few %, DECIGO 
and BBO will detect enough number of samples up to 
high redshifts. Including the lensing magnification, the 
Hubble parameter as a function of redshift up to z = 1 
can be determined with the accuracy of 3—14 %, 1.5—8 %, 
and 0.8 — 4%, for the observation time T b s = 1, 3, and 
10 yr, respectively. With a technique to remove the lens- 
ing magnification, we can further improve the accuracy 
by a factor of ~ 2.7 at z = 1. Although the feasibility 
of the method, particularly on the redshift indentifica- 
tion of the host galaxies and the detector calibration, 
still need to be investigated in more detail, the present 
method puts forward an interesting possibility to enlarge 
the science with standard sirens as a by-product, and 
it is complementary to other methodologies to measure 
Hubble parameter. 



8 



Acknowledgments 

We acknowledge N. Kanda, T. Tanaka, K. Yagi, J. 
Yokoyama, and C. Yoo for helpful discussions and valu- 
able comments. We also thank T. Smith for his effort 
to improve our manuscript. A. N. and S. S. are sup- 



ported by a Grant-in-Aid for JSPS Fellows. S. S. is also 
supported by JSPS through Excellent Young Researchers 
Overseas Visit Program. A. T. and S. S. is supported in 
part by a Grants-in-Aid for Scientific Research from the 
JSPS (No. 21740168 for A. T. and No. 21-00784 for S. 
S.). 



[1] N. Seto, S. Kawamura, and T. Nakamura, Phys. Rev. 

Lett. 87, 221103 (2001). 
[2] S. Sato et al., Journal of Physics: Conference Series 154, 

012040 (2009). 

[3] E. S. Phinney et al, The Big Bang Observer, NASA Mis- 
sion Concept Study (2003). 
[4] C. Cutler and D. E. Holz, Phys. Rev. D 80, 104009 

(2009) . 

[5] K. Yagi and T. Tanaka, Prog. Theor. Phys. 123, 1069 

(2010) . 

[6] A. Nishizawa, A. Taruya, and S. Kawamura, Phys. Rev. 

D 81, 104043 (2010). 
[7] B. F. Schutz, Nature (London) 323, 310 (1986). 
[8] C. Bonvin, R. Durrer, and M. Kunz, Phys. Rev. Lett. 

96, 191302 (2006). 
[9] C. Bonvin, R. Durrer, and A. Gasparini, Phys. Rev. D 

73, 023523 (2006). 
[10] D. Stern, R. Jimenez, L. Verde, M. Kamionkowski, and 

S. A. Stanford, J. Cosmol. Astropart. Phys. 02, 008 

(2010). 

[11] E. Gaztanaga, R. Miquel, and E. Sanchez, Phys. Rev. 

Lett. 103, 091302 (2009). 
[12] E. Gaztanaga, A. Cabre, and L. Hui, Mon. Not. R. As- 

tron. Soc. 399, 1663 (2009). 
[13] M. Sasaki, Mon. Not. R. Astron. Soc. 228, 653 (1987). 
[14] N. Jarosik et al., arXiv:1001.4744. 

[15] P. S. Corasaniti, M. L. Verde, A. Crotts, and C. Blake, 
Mon. Not. R. Astron. Soc. 369, 798 (2006). 

[16] A. Crotts et al, arXiv:astro-ph/0507043. 

[17] C. Cutler and E. E. Flanagan, Phys. Rev. D 49, 2658 
(1994). 

[18] M. Maggiore, Gravitational Waves (Oxford university 
press, 2008). 

[19] R. Takahashi and T. Nakamura, Prog. Theor. Phys. 113, 
63 (2005). 

[20] L. S. Finn, Phys. Rev. D 46, 5236 (1992). 
[21] T. A. Prince, M. Tinto, S. L. Larson, and J. W. Arm- 
strong, Phys. Rev. D 66, 122002 (2002). 
[22] C. Cutler and J. Harms, Phys. Rev. D 73, 042001 (2006). 
[23] R. Schneider, V. Ferrari, S. Matarrese, and S. F. P. 



Zwart, Mon. Not. R. Astron. Soc. 324, 797 (2001). 
[24] J. Abadie et al., Classical Quantum Gravity 27, 173001 
(2010). 

[25] J. Wambsganss, R. Cen, G. Xu, and J. P. Ostriker, As- 

trophys. J. 475, L81 (1997). 
[26] D. E. Holz and E. V. Linder, Astrophys. J. 631, 678 

(2005). 

[27] D. E. Holz and R. M. Wald, Phys. Rev. D 58, 063501 
(1998). 

[28] K. Kainulainen and V. Marra, arXiv:1011.0732. 

[29] C. M. Hirata, D. E. Holz, and C. Cutler, Phys. Rev. D 

81, 124046 (2010). 
[30] C. Shang and Z. Haiman, Mon. Not. R. Astron. Soc. 411, 

9 (2011). 

[31] C. Gordon, K. Land, and A. Slosar, Phys. Rev. Lett. 99, 
081301 (2007). 

[32] L. Silberman, A. Dekel, A. Eldar, and I. Zehavi, Astro- 
phys. J. 557, 102 (2001). 

[33] C. Shapiro, D. J. Bacon, M. Hendry, and B. Hoyle, Mon. 
Not. R. Astron. Soc. 404, 858 (2010). 

[34] S. Hilbert, J. R. Gair, and L. J. King, Mon. Not. R. 
Astron. Soc. 412, 1023 (2011). 

[35] Strictly speaking, this is true for flat geometry of the 
universe. For closed and open universe, H(z) is replaced 

with H{z) cos- 1 [/„' H^cosh- 1 [£ re- 

spectively. 

[36] Nevertheless, the gravitational lensing magnification 
might be statistically measurable and has been recently 
discussed in Ref. [|[ 

[37] For the geometric factor arising from the non-orthogonal 
detector arms (60°), we incorporate this effect into the 
detector noise curve P(f)- The factor arising from de- 
tector's angular response is also taken into account in 
evaluating the noise curve by averaging over the sky. 

[38] The noise power spectrum presented here is slightly dif- 
ferent from the one given by Ref. @, who considered or- 
thogonal detector arms and optimal incidence of a GW. 



